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Summary. Using the Rossi X-ray Timing Explorer (RossiXTE) , astronomers have 
discovered that disk-accreting neutron stars with weak magnetic fields produce three 
distinct types of high-frequency X-ray oscillations. These oscillations are powered 
by release of the binding energy of matter falling into the strong gravitational field 
of the star or by the sudden nuclear burning of matter that has accumulated in 
the outermost layers of the star. The frequencies of the oscillations reflect the or- 
bital frequencies of gas deep in the gravitational fleld of the star and/or the spin 
frequency of the star. These oscillations can therefore be used to explore fundamen- 
tal physics, such as strong-field gravity and the properties of matter under extreme 
conditions, and important astrophysical questions, such as the formation and evo- 
lution of millisecond pulsars. Observations using RossiXTE have shown that some 
two dozen neutron stars in low-mass X-ray binary systems have the spin rates and 
magnetic fields required to become millisecond radio-emitting pulsars when accre- 
tion ceases, but that few have spin rates above about 600 Hz. The properties of these 
stars show that the paucity of spin rates greater than 600 Hz is due in part to the 
magnetic braking component of the accretion torque and to the limited amount of 
angular momentum that can be accreted in such systems. Further study will show 
whether braking by gravitational radiation is also a factor. Analysis of the kilohertz 
oscillations has provided the first evidence for the existence of the innermost stable 
circular orbit around dense relativistic stars that is predicted by strong-field gen- 
eral relativity. It has also greatly narrowed the possible descriptions of ultradense 
matter. 



1 Introduction 



Neutron stars and black holes are important cosmic laboratories for studying 
fundamental questions in physics and astronomy, especially the properties of 
dense matter and strong gravitational fields. The neutron stars in low- mass 
X-ray binary systems (LMXBs) have proved to be particularly valuable sys- 
tems for investigating the innermost parts of accretion disks, gas dynamics 
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and radiation transport in strong radiation and gravitational fields, and the 

properties of dense matter, because the magnetic fields of many of these stars 
are relatively weak but not negligible. Magnetic fields of this size are weak 
enough to allow at least a fraction of the accreting gas to remain in orbit 
as it moves into the strong gravitational and radiation fields of these stars, 
but strong enough to produce anisotropic X-ray emission, potentially allow- 
ing the spin rates of these stars to be determined. The discovery and study 
using RossiXTE of oscillations in the X-ray emission of the accreting neutron 
stars in LMXBs, with frequencies comparable to the dynamical frequencies 
near these stars, has provided an important new tool for studying the strong 
gravitational and magnetic fields, spin rates, masses, and radii of these stars 
(for previous reviews, see [50, 51, 46, 96, 97]). 

Periodic accretion-powered X-ray oscillations have been detected at the 
spin frequencies of seven neutron stars with millisecond spin periods, estab- 
lishing that these stars have dynamically important magnetic fields. In this 
review, a pulsar is considered a millisecond pulsar (MSP) if its spin period Pg is 
< 10 ms (spin frequency i^spin > 100 Hz). The channeling of the accretion flow 
required to produce these oscillations implies that the stellar magnetic fields 
are greater than ~ 10^ G (see [69]), while the nearly sinusoidal waveforms of 
these oscillations and their relatively low amplitudes indicate that the fields 
are less than ~ 10^" G [78]. The spin frequencies of these accretion-powered 
MSPs range from 185 Hz to 598 Hz (see Table 1). 

Nearly periodic nuclear-powered X-ray oscillations (see Fig. 1) have been 
detected during the thermonuclear bursts of 17 accreting neutron stars in 
LMXBs, including 2 of the 7 known accretion-powered MSPs (Table 1). The 
existence of the thermonuclear bursts indicates that these stars' magnetic 
fields are less than ~ 10^" G [38, 54], while the spectra of the persistent X- 
ray emission [79] and the temporal properties of the burst oscillations (see 
[20, 18]) indicate field strengths greater than ^ lO"^ G. The spin frequencies 
of these nuclear-powered pulsars range from 45 Hz up to 1122 Hz. Three of 
them are also accretion-powered MSPs. 

Measurements of the frequencies, phases, and waveforms of the accretion- 
and nuclear-powered oscillations in SAX J1808.4— 3658 (see Fig. 1 and [20]) 
and XTE J1814— 338 [89] have shown that, except during the first seconds 
of some biirsts, the nuclear-powered oscillations have very nearly the same 
frequency, phase, and waveform as the accretion-powered oscillations, estab- 
lishing beyond any doubt (1) that these stars have magnetic fields strong 
enough to channel the accretion flow and enforce corotation of the gas at the 
surface of the star that has been heated by thermonuclear bursts and (2) that 
their nuclear- and accretion-powered X-ray oscillations are both produced by 
spin modulation of the X-ray flux from the stellar surface. The burst oscilla- 
tions of some other stars are also very stable [88], but many show frequency 
drifts and phase jitter [91, 92, 75, 73]. These results confirm that burst and 
persistent oscillations both reveal directly the spin frequency of the star. Sev- 
eral mechanisms for producing rotating emission patterns during X-ray bursts 
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Table 1. Accretion- and Nuclear-Powered Millisecond Pulsars" 





(Hz)'' 


Object 


Reference 


1122 


NK 


XTE J1739-285 


[431 


619 


NK 


4U 1608-52 


[291 


611 
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GS 1826-238 


[931 


601 


NK 


SAX J1750.8-2900 


[401 


598 


A 


IGR J00291+5934 


[611 


589 
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X 1743-29 


[861 


581 


NK 


4U 1636-53 


[109, 105, 92] 


567 
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MXB 1659-298 


1021 


549 


NK 


Aql X-1 


[1081 


530 
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A 1744-361 


[101 


524 


NK 


KS 1731-260 


[831 


435 
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XTE J1751-305 


[621 

L J 


410 


N 


SAX J1748.9-2021 


[41] 


401 


ANK 


SAX J1808.4-3658 


[103, 19] 


377 
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HETE J1900.1-2455 


[71] 


363 


NK 


4U 1728-34 


[91] 


330 


NK 


4U 1702-429 


[59] 


314 


AN 


XTE J1814-338 


[60] 


270 


N 


4U 1916-05 


[25] 


191 


AK 


XTE J1807.4-294 


[58, 101] 


185 


A 


XTE J0929-314 


[24] 


45 


N 


EXO 0748-676 


[99] 



"Defined in this review as pulsars with spin periods Ps < 10 ms. EXO 0748—676 is 
not a millisecond pulsar according to this definition. ''Spin frequency inferred from 
periodic or nearly periodic X-ray oscillations. A: accrotion-powcrcd millisecond 
pulsar. N: nuclear-powered millisecond pulsar. K: kilohcrtz QPO source. See text 
for details. Gavriil, Strohmayer, Swank & Markwardt have recently discovered a 
442 Hz X-ray pulsar in the direction of NGC 6440, but it is not yet certain whether 
it is nuclear- or accretion-powered. Consequently we have not counted this star in 
either category. 

have been proposed (see, e.g., [92, 72, 85, 8, 9, 26]), but which mechanisms 
are involved in which stars is not yet fully understood. 

Kilohcrtz quasi-periodic oscillations (QPOs) have now been detected in 
some two dozen accreting neutron stars (see [46]), including 9 of the 18 known 
nuclear-powered X-ray pulsars and 2 of the 7 known accretion-powered MSPs 
(Table 1). The frequencies of the kilohcrtz QPOs detected so far range from 
~1300 Hz in 4U 0614-^09 [98] down to ~10 Hz in Cir X-1 [17]. If, as expected, 
the frequencies of the highest-frequency kilohcrtz QPOs reflect the orbital 
frequencies of gas in the disk near the neutron star [69, 97], then in most 
kilohcrtz QPO sources gas is orbiting close to the surface of the star. The 
spin frequencies z^spin of the neutron stars in the kilohcrtz QPO sources are 
inferred from the periodic accretion- and nuclear-powered X-ray oscillations 
of these stars. In the systems in which kilohcrtz QPOs and periodic X-ray 
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Fig. 1. Left: X-ray burst and millisecond burst oscillations seen in 4U 1636—53. 
The main panel displays the X-ray countrate in a succession of 2-second time inter- 
vals, showing the rapid rise and approximately exponential decay of the burst. The 
inset panel shows the strong ~ 580 Hz X-ray countrate oscillations observed during 
the time interval bounded by the vertical dashed lines in the main panel. From T. 
Strohmayer, personal communication, see also [92]. Right: An X-ray burst with mil- 
lisecond X-ray countrate oscillations observed in SAX J1808.4— 3658 on 18 October 
2002. The dark curve and the scale at the right show the X-ray countrate as a func- 
tion of time during the burst. The contours show the dynamic power spectrum of 
the X-ray countrate on the scale at the left. Note the rapid increase in the oscillation 
frequency at the beginning of the burst, the disappearance of the oscillation at the 
peak of the burst, and its reappearance about 5 s later. The horizontal dashed line 
shows the frequency of the neutron star's spin inferred from its accretion-powered 
brightness oscillations. From [20]. 



oscillations have both been detected with high confidence (> 4ct), Jospin ranges 
from 191 Hz to 619 Hz. 

In many kilohertz QPO sources, the separation Z\i^qpo of the frequencies 
Uu and vg of the upper and lower kilohertz QPOs remains constant to within 
a few tens of Hz, even as and Vi vary by as much as a factor of 5. Avqpo 
is approximately equal to t'spin or i'spiii/2 in all stars in which these frequen- 
cies have been measured. In the accretion-powered MSPs XTE J1807.4— 294 
[55] and SAX J1808.4— 3658 [20], no variations of Z\i^QPo with time have so 
far been detected. In XTE J1807.4— 294, no difference between Z\i/qpo and 
I'spin has been detected; in SAX J1808.4— 3658, no difference between Z\i/qpo 
and i'spin/2 has been detected (see [55, 20, 96, 46, 104, 101]). These results 
demonstrate conclusively that at least some of the neutron stars that pro- 
duce kilohertz QPOs have dynamically important magnetic fields and that 
the spin of the star plays a central role in generating the kilohertz QPO pair. 
Consequently, Avqpo can be used to estimate, to within a factor of two, the 
otherwise unknown spin frequency of a star that produces a pair of kilohertz 
QPOs. 

The kilohertz QPO pairs recently discovered in Cir X-1 [17] extend sub- 
stantially the range of known kilohertz QPO behavior. In Cir X-1, values of I'u 
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Fig. 2. Power density spectrum of the variations in the X-ray countrate from the 
accretion-powered MSP SAX J1808.4-3658 seen on 18 October 2002. The peaks 
correspond to the 401 Hz periodic oscillations ("pulsations") at the star's spin fre- 
quency, the lower kilohcrtz QPO at 499±4 Hz, and the upper kilohcrtz QPO at 
694±4 Hz (from [104]). In this pulsar, the separation between the two kilohertz 
QPOs is half the spin frequency. Two kilohertz QPOs have also been seen in the 
accreting millisecond X-ray pulsar XTE J1807.4— 294, which has a spin frequency 
of 191 Hz [57, 101]. In this pulsar, the separation between the two kilohertz QPOs 
is consistent with the spin frequency. These results demonstrate that the star's spin 
plays a central role in the generation of kilohertz QPO pairs. 

and V(, as small, respectively, as 230 Hz and 50 Hz have been observed simul- 
taneously. These frequencies are more than 100 Hz lower than in any other 
kilohertz QPO system. Unlike the kilohertz QPO pairs so far observed in other 
neutron stars, in Cir X-1 Auq^o has been observed to increase with increasing 
Vu'- as Vu increased from ^230 Hz to ^500 Hz, Avqpq increased from 173 Hz to 
340 Hz. The relative frequency separations AvQPo/t^u in Cir X-1 are ~55%- 
75%, larger than the relative frequency separations ~20%-60% observed in 
other kilohertz QPO systems. Z\i^Qpo has been seen to vary by ^100 Hz in 
GX 5—1, which also has relatively low kilohertz QPO frequencies, but with no 
clear dependence on u^- If, as is generally thought, the frequencies of the kilo- 
hertz QPOs reflect the frequencies of orbits in the disk, kilohertz QPOs with 
such low frequencies would require the involvement of orbits ~50 km from the 
star, which is a challenge for existing kilohertz QPO models. Accretion- and 
nuclear-powered X-ray oscillations have not yet been detected in Cir X-1 and 
hence its spin frequency has not yet been measured directly. Further study of 
Cir X-1 and the relatively extreme properties of its kilohertz QPOs is likely 
to advance our understanding of the physical mechanisms that generate the 
kilohertz QPOs in all systems. 

The first 11 spins of accretion-powered X-ray MSPs that were measured 
were consistent with a flat distribution that ends at 760 Hz [20] , but were also 
consistent with a distribution that decreases more gradually with increasing 
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frequency [66] ; the spins of 22 accretion- and nuclear-powered X-ray MSPs are 
now known. The proportion of accretion- and nuclcar-powcrcd MSPs with fre- 
quencies higher than 500 Hz is greater than the proportion of known rotation- 
powered MSPs with such high frequencies, probably because there is no bias 
against detecting X-ray MSPs with high frequencies, whereas detection of 
rotation-powered radio MSPs with high spin frequencies is still difHcult [20]. 
The recent discovery of a 1122 Hz MSP [43] supports this argument, which is 
not inconsistent with the recent discovery of a 716 Hz rotation-powered radio 
MSP [31]. 

These discoveries have established that many neutron stars in LMXBs have 
magnetic fields and spin rates similar to those of the rotation-powered MSPs. 
The similarity of these neutron stars to rotation-powered MSPs strongly sup- 
ports the hypothesis [2, 80] that they are the progenitors of the rotation- 
powered MSPs. After being spun down by rotation-powered emission, the 
neutron stars in these systems are thought to be spun up to millisecond pe- 
riods by accretion of matter from their binary companions, eventually be- 
coming nuclear- and accretion-powered MSPs and then, when accretion ends, 
rotation-powered MSPs. 

In §2 we discuss in more detail the production of accretion- and rotation- 
powered MSPs by spin-up of accreting weak-field neutron stars in LMXBs, 
following [52], and in §3 we describe several mechanisms that may explain the 
nuclear-powered X-ray oscillations produced at the stellar spin frequency by 
such stars. In §§4 and 5 we discuss, respectively, possible mechanisms for gen- 
erating the kilohertz QPO pairs, following [49], and how the kilohertz QPOs 
can be used as tools to explore dense matter and strong gravity, following 
[69, 50, 51]. 

2 Production of Millisecond Pulsars 

Neutron stars in LMXBs are accreting gas from a Keplerian disk fed by a low- 
mass companion star. The star's magnetic field and accretion rate are thought 
to be the most important factors that determine the accretion flow pattern 
near it and the spectral and temporal characteristics of its X-ray emission (see 
[69]). The accretion rates of these stars vary with time and can range from the 
Eddington critical rate Me to less than 10~^M£;. Their magnetic flelds are 
thought to range from 10^^ G down to 10^ G or possibly less, based on their 
X-ray spectra [79], the occurrence of thermonuclear X-ray bursts [38], and 
their high-frequency X-ray variability [69, 96]. Magnetic fields at the upper 
end of this range are strong enough to terminate the Keplerian disk well above 
the stellar surface, even for accretion rates ^Me, whereas magnetic fields at 
the lower end of this range affect the flow only close to the star, even for 
accretion rates as low as ~10~*M£;. 

For intermediate field strengths and accretion rates, some of the accreting 
gas is expected to couple to the star's magnetic field well above the stellar 
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Fig. 3. Side view of a weak-field neutron star accreting from a disk, showing the 
complex flow pattern expected. Some accreting gas couples strongly to the magnetic 
field and is channeled toward the magnetic poles, but a substantial fraction couples 
only weakly and drifts inward in nearly circular orbits as it transfers its angular 
momentum to the star via the stellar magnetic fleld. From [69]. 



surface. The star's magnetic field channels part of the flow toward its magnetic 
poles, and this flow heats the outer layers of the star unevenly. The remainder 
of the accreting gas is expected to remain in a geometrically thin Keplerian 
flow that penetrates close to the stellar surface, as shown in Figure 3. The 
gas that remains in orbit close to the star is thought to be responsible for 
generating the kilohertz QPOs (see [48, 46, 49]). When thermonuclear X-ray 
bursts occur, they also heat the outer layers of the star unevenly. Whether 
due to accretion or to nuclear burning, the uneven heating of the outer layers 
produces a broad pattern of X-ray emission that rotates with the star, making 
both the accretion-powered and nuclear-powered X-ray emission of the star 
appear to oscillate at the spin frequency. The stability of the nuclear-powered 
oscillations show that the heated region is strongly coupled to the rotation 
of the star, probably via the star's magnetic field. The phase locking of the 
nuclear- and accretion-powered oscillations and the strong similarity of the 
two waveforms in SAX J1808.4-3658 and XTE J1814-338 indicate that the 
stellar magnetic field is playing a dominant role, at least in these pulsars. 
However, these two are the only nuclear-powered pulsars in which accretion- 
powered oscillations at the spin frequency have also been detected. 

Production of millisecond accretion-powered pulsars. — Accretion from a 
disk will spin up a slowly-rotating neutron star on the spin-relaxation timescale 
[27, 28, 52] 

t . = 27r^,pinJ / z/,pi„ ^^ M \ 

^■'^ - [M(GA/r™)V2] IsOO Uz) [o.OIMe ) ' ^ ^ 

where Jospin, M , and / are the star's spin rate, mass, and moment of inertia, 
M is the accretion rate onto the star (not the niRSS trRiisfcr rate), r,^^ is the 
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angular momentum coupling radius, a is 0.23 if the inner disk is radiation- 

prcssurc-dominatcd (RPD) or 0.38 if it is gas-prcssurc-dominated (GPD), and 
in the last expression on the right the weak dependence of tg^m on M, /, and 
the star's magnetic field has been neglected. 

The current spin rates of neutron stars in LMXBs depend on the average 
accretion torque acting on them over a time ~ tspin- Determining this average 
torque is complicated by the fact that the accretion rates and magnetic fields 
of these stars vary with time by large factors and that the accretion torque 
can decrease as well as increase the spin rate. Mass transfer in the neutron- 
star-white-dwarf binary systems is thought to be stable, with a rate that 
diminishes secularly with time. 

While a few neutron stars in LMXBs accrete steadily at rates ~ Me-, most 
accrete at rates ~10~^-10~^M£; or even less [30, 45, 95, 69] and many ac- 
crete only episodically [95, 81]. Important examples are the known accretion- 
powered MSPs in LMXBs, which have outbursts every few years during which 
their accretion rates rise to ~10~^M£; for a few weeks before falling again to 
less than ^X^^^Me [20, 89]. Also, there is strong evidence that the external 
magnetic fields of neutron stars in LMXBs decrease by factors ~10^-10^ dur- 
ing their accretion phase, perhaps on timescales as short as hundreds of years 
(see [82, 7]). 

If a star's magnetic field and accretion rate are constant and no other 
torques are important, accretion will spin it up on a timescale ~ fspin to its 
equilibrium spin frequency t'eq. This frequency depends on A'/, the strength 
and structure of the star's magnetic field, the thermal structure of the disk at 
Tm, and M [27, 100, 28]. If a star's magnetic field and accretion rate change on 
timescales longer than tgpin, the spin frequency will approach t'eq and track it 
as it changes. If instead M varies on timescales shorter than igpin, the spin rate 
will fluctuate about the appropriately time-averaged value of feq (sec [23]). 
Thus feq and its dependence on B and M provide a framework for analyzing 
the evolution of the spins and magnetic fields of neutron stars in LMXBs. 

Figure 4 shows i^oq for five accretion rates and dipolc magnetic fields B^, 
assumed given by 3.2 X 10i9(PP)i/2 G and ranging from 10^ G to lO^'^ G. 
The lines are actually bands, due to systematic uncertainties in the models. 
The lines for M = Me and M = O.IM^; have jumps where the structure of 
the disk at the angular momentum coupling radius changes from RPD 
(lower left) to GPD (upper right); in reality the transition is smooth. For M 
less than ^ O.OIM^, the disk is GPD at even if the star's magnetic field 
is less than ~ 3 x 10^ G. Not shown are the effects of the stellar surface and 
the innermost stable circular orbit [52], which affect the spin evolution at spin 
periods less than ~ 1 ms. 

The spin rates of the known MSPs in LMXBs (see Table 1) are consistent 
with spin-up by accretion. The existence of only a single candidate with a 
spin rate greater than 620 Hz could be because (1) these stars have reached 
accretion spin equilibrium and fgq is less than 620 Hz for their appropriately 
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Fig. 4. Spin-evolution diagram. Lines sloping downward to the right show the P-P 
relation for magnetic dipole braking by a field with the strength indicated. Lines 
sloping upward to the right show the equilibrium spin period of a neutron star with 
the accretion rate indicated by the labels and a dipole field of the strength indicated 
by the downward-sloping lines. The dashed line sloping upward to the right shows 
where stars with a spin-down time equal to 15 Gy would lie. Data points are known 
rotation-powered pulsars; those of pulsars in binary systems are encircled. From [52]; 
data from [35]. 

(torque- weighted) time-averaged accretion rates, (2) they are still spinning up 
but the spin-up timescales for their current accretion rates are longer than the 
times they have been accreting at these rates, or (3) an additional braking 
torque is acting on them. 

For example, the 45 Hz spin rate of the pulsar EXO 0748—676 corresponds 
to accretion spin equilibrium for a dipole magnetic field of 2 x 10^ G and a 
time-averaged accretion rate of 10~^ Me, giving a spin-evolution time scale of 
~20 Myr, whereas the 191 Hz spin rate of XTE J1807. 4-294 corresponds to 
equilibrium for a field of 3 x 10^ G and M « 10~'^Me, giving a spin-evolution 
time scale '^80 Myr. For comparison, the 600 Hz spin rate of 4U 1608—52 cor- 
responds to equilibrium for a field of 3 x 10'' G and a time-averaged accretion 
rate of 10~^ Me, giving a spin-evolution time scale ^2 Gyr. These examples 
show that the spin rates of the known MSPs in LMXBs are consistent with 
spin-up to accretion spin equilibrium if they have magnetic fields in the range 
^3 x 10'' G to ~2 X 10^ G and average accretion rates in the range ^10~^ Me 
to ~10"^M£;, but that stars with accretion rates less than ^10~^Me may 
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not be in spin equilibrium but instead spinning up on a timescale longer than 

their accretion phases. In particular, the number of MSPs with spin frequen- 
cies > 620 Hz [20] may be small because the equilibrium spin rates of these 
stars are < 620 Hz or because their spin-up timescales are longer than their 
accretion phases. As an example, the timescale to reach the 1122 Hz spin fre- 
quency reported for XTE J1739— 285 [43] is 400 Myr for a long-term average 
accretion rate of 10~^Me and 3Gyr for an accretion rate of 10~^Me. The 
ranges of magnetic fields and accretion rates required are consistent with the 
other observed properties of neutron stars in LMXBs [69, 78, 20]. 

If their magnetic fields are weak enough, it is possible that the spin rates of 
some neutron stars in LMXBs are affected by gravitational radiation torques. 
Based on the limited information then available, some authors [13, 94] spec- 
ulated that neutron stars in LMXBs have negligible magnetic fields and spin 
frequencies in a narrow range, with many within 20% of 300 Hz. Such a dis- 
tribution would be difficult to explain by accretion torques and was taken as 
evidence that gravitational radiation plays an important role. We now know 
(sec § 2) that most if not all neutron stars in LMXBs have dynamically im- 
portant magnetic fields, that the observed spins of neutron stars in LMXBs 
are distributed roughly uniformly from < 200 Hz to > 600 Hz, and that pro- 
duction of gravitational radiation by uneven heating of the crust or excitation 
of r- waves is not as easy as was originally thought [94, 56] . At present there is 
no unambiguous evidence that the spin rates of neutron stars in LMXBs are 
affected by gravitational radiation. 

Production of millisecond rotation-powered pulsars. — Soon after rotation- 
powered radio-emitting MSPs were discovered, it was proposed that they have 
been spun up to millisecond periods by steady accretion in LMXBs at rates 
~ Me (see [12]), with the implicit assumption that accretion then ends sud- 
denly; otherwise the stars would track feq to low spin rates as the accretion 
phase ends. This simplified picture is sometimes still used (sec, e.g., [3]), but — 
as noted above — most neutron stars in LMXBs accrete at rates <C Me, many 
accrete only episodically, and the accretion rates of others dwindle as their 
binary systems evolve. The real situation is therefore more complex. 

The initial spins of rotation-powered MSPs recycled in LMXBs are the 
spins of their progenitors when they stopped accreting. These spins depend 
sensitively on the magnetic fields and the appropriately averaged accretion 
rates of the progenitors when accretion ends. Comparison of the equilibrium 
spin-period curves for a range of accretion rates with the P-P distribution of 
known rotation-powered MSPs (Fig. 4) suggests three important conclusions: 

(1) The hypothesis that the accretion torque vanishes at a spin frequency 
close to the calculated i^cq predicts that MSPs should not be found above the 
spin-equilibrium line for M = Me, because this is a bounding case. The ob- 
served P-P distribution is consistent with this requirement for the RPD model 
of the inner disk that was used for M greater that ~ O.IM^, except for two 
pulsars recently discovered in globular clusters: B1821— 24 and B1820— 30A 
[34]. Either the intrinsic P's of these pulsars are lower than shown in Fig. 4, 
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or the RPD model of the inner disk does not accurately describe the accretion 
flow that spun up these stars. 

(2) The accretion spin-equilibrium hypothesis predicts that MSPs should 
be rare or absent below the spin-equilibrium line for M = 10~'^Me, because 
stars accreting at such low rates generally will not achieve millisecond spin 
periods during their accretion phase. The observed P-P distribution is con- 
sistent with this prediction. 

(3) The MSPs near the 15 Gyr spin-down line were produced in situ by 
final accretion rates less than ~ 3 x 10~^Me rather than by spin-up to shorter 
periods by accretion at rates greater than ~ 3 x 10~^Me followed by magnetic 
braking, because braking would take too long. This result accords with the 
expectation (see above) that most neutron stars in LMXBs accrete at rates 
<C Me toward the end of their accretion phase. 

3 Nuclear-Powered X-ray Oscillations 

Accretion of matter onto the surface of a neutron star produces a fluid ocean 
on top of the solid crust (see [84]). Depending on the accretion rate and the 
initial composition of the material, the conditions needed to ignite hydrogen 
and helium can be reached. Ignition of the accreted matter will generally occur 
at a particular place on the surface of the star. For low to moderate accretion 
rates, burning is unstable and produces type-I X-ray bursts (see [85]). There 
are several important timescales in this problem, including the time required 
for burning to spread over the surface of the star, the time required for heat 
to reach the photosphere, and the timescale on which the heated matter cools 
by emission of radiation. The time required for burning to spread is expected 
to depend on the latitude (s) of the ignition point (s), because of the variation 
with latitude of the Coriolis force, which affects the thickness of the burning 
front and hence the speed at which it advances. A recent simulation [84] finds 
that the spreading time is shorter if the ignition point is nearer the equator, 
because there the burning front is less steep and propagates faster. The time 
required for burning to spread around the star is expected to be less than a 
second [11], much smaller than the ^ 10-30 s observed durations of the bursts 
in X-rays, which may reflect the time required for heat from the burning layers 
to reach the photosphere. If so, nuclear burning is probably over, or almost 
over, by the time the burst becomes visible in X-rays to a distant observer. 

Useful information about the burst evolution can be obtained from the 
nearly coherent X-ray oscillations seen during portions of some bursts (see 
[85]). The discovery of burst oscillations with very nearly the same frequen- 
cies as the spin rates of two MSPs [20, 89], as well as the observed stability of 
these oscillation frequencies [88] , eliminated any doubt that burst oscillations 
are generated by the spin of the star. However, several important questions 
are not yet fully resolved. In most bursters, the oscillation frequencies vary 
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slightly, especially during the burst rise, but in burst tails the oscillation fre- 
quencies often approach an asymptotic value [87] that remains the same to 
high precision over long times for a given source [90] . Determining what pro- 
duces these oscillations and what causes the differences in their behavior from 
burst to burst and star to star is important for understanding the physics of 
the bursts. 

The most widely discussed picture for type-I bursts assumes a hotter region 
on the surface of the star that has been heated from below and rotates with 
the star. The increase in the oscillation frequency observed near the beginning 
of some bursts has been attributed to the collapse of the stellar atmosphere 
that would occur as it cools after the end of nuclear burning [86, 21]. In this 
model, the rotation rate of the outer envelope and photosphere increases as 
the outer layers collapse at approximately constant angular momentum. This 
model is believed to capture an important aspect of the actual rotational 
behavior of the envelope during the rise of a burst, even though the observed 
frequency changes are larger than those predicted by the model by factors 
2- 3 [22] and it is not clear how uniform rotation of the gas in the envelope 
can be maintained in the presence of Coriolis, magnetic, and viscous forces. 
During the ~0.1 s rise of bursts, oscillations with relative amplitudes as high 
as 75% are observed [92] with frequencies that differ from the stellar spin 
frequency by up to ~ 1% [25]. This model probably does not provide a good 
desciption of the oscillations during the tails of bursts, when the temperature 
of the stellar surface is expected to be relatively uniform. During the burst 
tails, stable X-ray oscillations with amplitudes as large as 15% are observed 
for up to 30 s, with frequencies that are consistent, within the errors, with the 
stellar spin frequency [76, 20, 89]. The amplitudes of the oscillations during 
the rise of bursts appear to be anticorrelated with the X-ray flux, whereas 
no such relation is apparent during the tails of bursts [76] . These differences 
suggest that different mechanisms are responsible for the oscillations near the 
beginning and in the tails of bursts. 

Excitation by bursts of r- waves and other nonradial modes (see, e.g., [63]) 
in the surface layers of neutron stars has been proposed [32, 33] as a possible 
mechanism for producing observable X-ray oscillations during bursts. This 
idea has been explored further [53, 77], still without distinguishing between 
the oscillations observed near the beginnings of bursts and those observed 
in the tails of bursts. As noted above, these have significantly different char- 
acteristics. An important challenge for oscillation mode models is to explain 
the relatively large frequency variations observed during bursts. It has been 
suggested that these variations can be explained by changes in the character 
of the oscillations during the burst (e.g., from an r-mode to an interface mode 
[32]). Other challenges for models that invoke oscillations in the surface layers 
of the neutron star are to explain what physics singles out one or a few modes 
from among the very large number that could in principle be excited and how 
these modes can produce the X-ray oscillations with the relatively large am- 
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Fig. 5. Left: Correlation between the upper and lower kilohertz QPO frequencies 
seen in all sources in which both have bcou detected. The filled stars indicate fre- 
quencies seen in Atoll sources; the empty circles indicate those seen in Z sources. 
The shaded boxes (bottom left) indicate the frequencies of the kilohertz QPO pair 
recently discovered in Cir X-1. The shaded stars indicate the frequencies of the QPO 
pair seen in the MSP SAX J1808.4-3658. A few high frequency QPOs from black 
hole systems are included; the various lines represent fits of formulas to parts of the 
data. Prom [15]. Right: Quality factors of the lower and upper kilohertz QPOs ob- 
served in 4U 1636—53 and the expected frequency-dependence of the quality factor 
predicted by a model of an active oscillating region approaching the ISCO. Prom [4]. 



plitudes and high coherence observed. Further work is needed to resolve these 
questions. 

4 Accretion-Powered Kilohertz QPOs 

The properties of the kilohertz QPO pairs provide strong hints about the 
mechanisms that generate them (for a more complete discussion, see [49]): 

1. It appears very likely that the frequency of one of the two kilohertz 
QPOs reflects the orbital frequency of gas in the inner disk. The frequencies of 
the kilohertz QPOs are similar to those of orbital motion near neutron stars. 
They also vary by hundreds of Hertz on time scales as short as minutes (see, 
e.g., [65, 96, 101]). Such large, rapid variations are possible if they are related 
to orbital motion at a radius that varies [46]. 

2. The star's spin is somehow involved in producing the frequency sep- 
aration of the two kilohertz QPOs in a pair. This involvement is clear in 
XTE J1807-294, where Z\i^QPo ~ i^spin, and in SAX J1808.4-3658, where 
^J^QPO ~ i'spin/2. It is strongly indicated in the other kilohertz QPO sources, 
because in all cases where both Auqpq and fgpin have been measured, the 
largest value of Z\i^qpo is consistent or approximately consistent with either 

^spin 

or j/,pin/2 (see [96, 46, 49]). 



14 EYederick K. Lamb and Stratos Boutloukos 

3. A mechanism that produces a single sideband is indicated. Most mech- 
anisms that modulate the X-ray brightness at two frequencies (siich as am- 
plitude modulation) would generate at least two strong sidebands. Although 
weak sidebands have been detected close to the frequency of the lower kilo- 
hertz QPO [36, 37], at most two strong kilohcrtz QPOs are observed in a given 
system [96, 64]. This suggests that the frequency of one QPO is the primary 
frequency while the other is generated by a single- sideband mechanism. Beat- 
frequency mechanisms naturally produce a single sideband. Because one QPO 
frequency is almost certainly an orbital frequency, the most natural mecha- 
nism would be one in which the second frequency is generated by a beat with 
the star's spin frequency. 

4. Mechanisms for generating kilohertz QPO pairs like the 3:2 resonance 
proposed to explain the high-frequency QPOs observed in black hole can- 
didates [1] and the kilohertz QPOs observed in SAX J1808.4-3658 [44] are 
excluded as explanations for the kilohertz QPO pairs seen in neutron stars, 
because these mechanisms require a low-order resonance between the geodesic 
frequencies of test particles orbiting at a fixed radius, which disappears when 
the two frequencies change substantially, as they do in neutron stars (see, e.g., 
Fig. 5). As noted above, in many neutron stars the separation frequency is 
approximately constant, which is incompatible with a fixed frequency ratio 
[14, 15]. This type of mechanism also cannot explain the commensurability 
of Auqpo with z/gpin in all neutron stars in which both frequencies have been 
measured [47]. 

5. Production of kilohertz QPOs by oscillating emission from a narrow 
annulus in the accretion disk is incompatible with their quality factors and 
amplitudes, for the following reason [49]: A kilohertz QPO peak of relative 
width ^i^Qpo/i^QPO corresponds to the spread of geodesic frequencies in an 
annulus of relative width dr/r ~ ^i'qpo/'^qpo- The emission from an annulus 
in the inner disk of relative width Sr/r is a fraction ^ (5r/r)Ldisk of the 
emission from the entire disk and hence a fraction ~ (5r/r) [.i^disk/(-Ldisk + 
Lstur)] of the emission from the system. Thus the relative amplitude of a 
QPO of width SuQpo produced by oscillating emission from such an annulus 
is ^ ('^I'QPo/i'QPo) [^disk/(idisk + istar)]- Some kilohertz QPOs have relative 
widths Svqpo/vqpo % 0.005 (see [6, 69, 96, 97, 4, 5, 101]) and the accretion 
luminosity of a neutron star is typically ^ 5 times the accretion luminosity of 
the entire disk [67] . Consequently, even if the emission from the annulus were 
100% modulated at the frequency of the kilohertz QPO, which is very unlikely, 
the relative amplitude of the QPO would be only ^ 0.005 x 1/6 ^ 0.08%, much 
less that the 2-60 keV relative amplitudes ~ 15% observed in many kilohertz 
QPO sources (see, e.g., [69, 96, 97, 101]). 

A recently proposed modification [49] of the original sonic-point beat- 
frequency model [69] potentially can explain within a single framework why 
the frequency separation is close to z/gpin in some stars but close to i'spin/2 
in others. In this "sonic-point and spin-resonance" (SPSR) beat-frequency 
model, gas from perturbations orbiting at the sonic-point radius Vgp produces 
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a radiation pattern rotating with a frequency i^u close to the orbital frequency 
forb at ^spi SiS in the original model, and this rotating pattern is detected as 
the upper kilohertz QPO. This mechanism for generating the upper kilohertz 
QPO is supported by the observed anticorrelation of the upper kilohertz QPO 
frequency with the normal branch oscillation flux in Sco X-1 [107] and the 
anticorrelation of the kilohertz QPO frequency with the mHz QPO flux in 
4U 1608-52 [106]. 

A new ingredient in the modiflcd model is preferential excitation by the 
magnetic and radiation fields rotating with the neutron star of vertical motions 
in the disk at the "spin-resonance" radius rgr where i^spin — t'orb is equal to the 
vertical epicyclic frequency ly-^. Preliminary numerical simulations show that 
the resulting vertical displacement of the gas in the disk is much greater at the 
resonant radius than at any other radius. In a Newtonian 1/r gravitational 
potential, v^{r) — i'orb(^)- Although f^pir) is not exactly equal to Vo-chi''') 
in general relativity, the difference is < 2 Hz at the radii of interest (where 
^orh < 300 Hz). Consequently, at the resonance radius where vertical motion 
is preferentially excited, I'orb ~ z^^ ~ i^spin/2. At this radius, the orbital and 
vertical frequencies are both approximately i'spin/2. 

In the SPSR model, the clumps of gas orbiting the star at the sonic radius 
Tsp act as a screen, forming the radiation from the stellar surface into a pat- 
tern that rotates around the star with frequency forb(?'sp)- Interaction of this 
rotating radiation pattern with the gas in the disk that has been excited ver- 
tically at Tsr produces a second QPO with frequency vi = Vorh{fsp) — t'spin/2, 
if the gas at rgr is highly clumped, or with frequency = Vorhifsp) " t'spim 
if the flow at Vsr is relatively smooth. This second QPO is the lower kilohertz 
QPO. 

To see how the observed QPO frequency relations can be generated, sup- 
pose flrst that the distribution of the gas in the disk near the spin-resonance 
radius is relatively smooth. There may be a large number of small clumps 
or the flow may even be smooth. Each element of gas is oscillating vertically 
with frequency i'spin/2, but together they form a pattern of raised fluid el- 
ements that rotates around the star with frequency Vs\yn\- Becaiisc a large 
number of fluid elements are scattering radiation to the observer at any given 
moment, their individual contributions blend together, so the dominant time 
variation has frequency forb(^sp) — ^'spin- In this case the brightness variation 
produced by the pattern of scattering clumps dominates the brightness vari- 
ation produced by the individual clumps. The left-hand panel of Fig. 6 shows 
the power spectrum of the flux variation generated in a simulation in which 
500 randomly-positioned clumps scatter the radiation pattern coming from 
the sonic radius. The peak at Uorhii'sp) — Jospin is clearly dominant. 

Suppose instead that the gas in the disk near the spin-resonance radius is 
highly clumped. When illuminated, each clump orbiting at Vsr scatters radia- 
tion in all directions. In eff'ect, each clump redirects the radiation propagating 
outward from the sonic radius in the modest solid angle that it subtends (as 
seen from the sonic radius) into all directions. From the point of view of a 
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Fig. 6. Power spectra of the X-ray flux modulation produced by simulations of 
a disk with a large number of clumps near the spin-resonance radius (left-hand 

panel) and a small number of clumps (right-hand panel). The star's spin frequency 
is 400 Hz while the orbital frequency at the sonic radius is 1,000 Hz. These power 
spectra demonstrate that if the flow near the spin-resonance radius is relatively 

smooth, the effect of the clump pattern dominates and the dominant frequency 
is i^orb('^'sp) — i^spin. If iiistcad the flow is clumpy, the effect of individual clumps 



dominates and the dominant frequency is forb(rsp) 



[i/2. This simulation did 



not include any signal with the orbital frequency of the gas at the sonic radius. 



distant observer, each individual clump looks like a light bulb that is blinking 
on and off with a frequency equal to /^rb(fsp) — '^orb('fsr) ~ '^rb(fsp) — '^spin/Z. 
If there are only a modest number of clumps at rsr, the scattering from the 
individual clumps dominates the time variation of the X-ray flux. The right- 
hand panel of Fig. 6 shows the power spectrum of the flux variation generated 
in a simulation in which five randomly-positioned clumps scatter the radia- 
tion pattern coming from the sonic radius. The peak at i'orb(?'sp) — i^spin/S 
is clearly dominant. Because the radiation is scattered in all directions, an 
observer does not have to be close to the disk plane to see the X-ray flux 
modulation. 

Magnetic forces may cause the gas in the accretion disk to become more 
clumped as it approaches the neutron star [69, 48, 49]. Consequently, the 
parameters that may be most important in determining whether the flow at 
the spin resonance radius rsr is clumpy or smooth are the star's spin frequency 
and magnetic field. For a given stellar magnetic field, the flow is likely to be 
more clumpy if the star is spinning rapidly and rgr is therefore close to the 
star. For a given spin rate, the flow is likely to be more clumpy if the star's 
magnetic fleld is stronger. 

The four sources with i^spin < 400 Hz and measurable frequency sepa- 
rations have AfQpo « Vspin whereas the five sources with i/gpin > 400 Hz 
have Ai/Qpo ~ t'spin/2 (see [74]). With such a small sample, one cannot make 
any definite statements, but the apparent trend is consistent with the sonic- 



Accreting Neutron Stars in Low-Mass X-Ray Binary Systems 



17 



point and spin-resonance beat-frequency model. These trends suggest that if 
kilohcrtz QPOs are detected in the recently-discovered 185 Hz and 314 Hz 
accretion-powered X-ray pulsars XTE J0929-314 [24] and XTE J1814-338 
[89], their frequency separations should be approximately equal to their re- 
spective spin frequencies. The 435 Hz spin frcriucncy of XTE .71751—305 [62] 
is high enough that Avqpq could be either approximately 435 Hz or approx- 
imately 217 Hz; QPOs at both frequencies might even be detectable. 

Finally, we note that there is no known reason why the mechanism for 
producing a lower kilohertz QPO proposed in the original sonic-point beat- 
frequency model would not operate. Apparently this mechanism does not 
produce a strong QPO in the fast rotators, but it might produce a weak QPO 
in these sources. If it operates in the slow rotators, it would produce a QPO 
near forb(''sp) — J'spin that might appear as a sideband to the lower kilohertz 
QPO. 

The sonic-point and spin-resonance beat-frequency model appears qualita- 
tively consistent with the basic properties of the kilohertz QPOs, but whether 

it can explain their detailed properties and the wide range of frequencies seen 
in different systems, such as Circinus X-1, remains to be determined. 

5 Kilohertz QPOs as Tools 

As explained in the previous section, despite uncertainty about the precise 
physical mechanisms responsible for generating the kilohertz QPO pairs seen 
in neutron star systems, there is good evidence that the upper kilohertz QPO 
is produced by orbital motion of gas in the strong gravitational field near 
the star. Making only this minimal assumption, the kilohertz QPOs can be 
used as tools to obtain important constraints on the masses and radii of the 
neutron stars in LMXBs and explore the properties of ultradense matter and 
strong gravitational fields (see [69, 50, 51]). 

For example, the left panel of Fig. 7 shows how to construct constraints 
on the mass and radius of a nonrotating neutron star, given i/*, the highest 
orbital frequency observed in the source. i?orb must be greater than the stellar 
radius, so the star's representative point must lie to the left of the (dashed) 
cubic curve M°(i?orb) that relates the star's mass to the radius of orbits 
with frequency f*. The high coherence of the oscillations constrains i?orb to 
be greater than Rms, the radius of the innermost stable orbit, which means 
that the radius of the actual orbit must lie on the M°(i?,orb) curve below its 
intersection with the (dotted) straight line M°(i?ins) that relates the star's 
mass to Rms- These requirements constrain the star's representative point to 
lie in the unhatched, pie-slice shaped region enclosed by the solid line. The 
allowed region shown is for v* = 1330 Hz, the highest value of i/„ observed in 
4U 0614-1-09 [98], which is also the highest value so far observed in any source. 

The right panel of Fig. 7 shows how this allowed region c;oniparc!s with 
the mass-radius relations given by five representative equations of state (for 
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Fig. 7. Radius-mass plane showing the constraints on neutron star masses and 

radii and the equation of state of neutron-star matter that can be derived from the 
frequency of the upper kilohertz QPO, here 1330 Hz, which is thought to be the 
orbital frequency of gas accreting onto the star. Left panel: The dashed curved line 
shows the relation between the mass of the star and the radius -Rorb of the orbit for 
a nonrotating star, which is an upper bound on the radius of a nonrotating star. The 
diagonal dotted line shows the relation between the mass of the star and the radius 
J?ms of the marginally stable orbit, which must be larger than Roih in order for the 
gas to make the hundreds of orbits around the star indicated by the coherence of 
the kilohertz QPO waveform. Consequently the mass and radius of the star must 
correspond to a point inside the unshaded "slice of pie". If the QPO frequency is 
shown to be that of the marginally stable orbit, then Rorb = Rms and the mass 
of the star is determined precisely. Right panel: Curves of mass-radius relations for 
nonrotating stars constructed using several proposed neutron-star matter equations 
of state (EOS), showing that a 1330 Hz QPO is just inconsistent with EOS M. The 
higher the observed QPO frequency, the tighter the constraints. After [69]. 

a description of these EOS and references to the literature, see [68]). If 
4U 0614+09 were not spinning, EOS L and M would both be excluded. How- 
ever, 4U 0614+09 is undoubtedly spinning (the frequency separation Avqpo 
between its two kilohertz QPOs varies from 240 Hz to 355 Hz [98]). If its spin 
frequency is high, EOS M may be allowed, but EOS L is excluded for any spin 
rate. 

Assuming that the upper kilohertz QPO at Uu is produced by orbital mo- 
tion of gas near the neutron star, its behavior can be used to investigate orbits 
in the region of strongly curved spacetime near the star. For example, it may 
be possible to establish the existence of an innermost stable circular orbit 
(ISCO) around some neutron stars in LMXBs (see [39, 69, 51]. This would 
be an important step forward in determining the properties of strong gravi- 
tational fields and dense matter, because it would be the first confirmation of 
a prediction of general relativity in the strong-field regime. 
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The sonic-point model of the kilohertz QPOs predicts several signatures 
of the ISCO [68, 50]. As an example, it predicts that the frequencies of both 
kilohertz QPOs will increase with the increasing accretion luminosity until 
the sonic radius — which moves inward as the mass flux through the inner disk 
increases — reaches the ISCO, at which point the frequencies of the kilohertz 
QPOs will become approximately independent of the accretion luminosity. 
Behavior similar to this has been observed [110, 42, 16], but important issues, 
such as the robustness of the predicted relation between QPO frequency and 
M, need further work. 

The sonic-point model also predicts a steep drop in the coherence of the 
kilohertz QPOs as the orbits involved approach the ISCO [50, 69, 70]. Abrupt 
drops have been observed in the quality factors of the kilohertz QPOs in 
several atoll sources, consistent with models of their expected behavior as the 
orbit involved approaches the ISCO [4, 5]. 

If either of these behaviors can be shown to be caused by an ISCO, it will 
be a major advance in establishing the properties of strong-field gravity. 

We thank D. Chakrabarty, C.J. Cook, J.M. Cook, M. van dcr Khs, M.C. 
Miller, and J. Swank for helpful discussions. This research was supported in 
part by NASA grant NAG5-12030, NSF grant AST0098399, and funds of the 
Fortner Endowed Chair at Illinois. 
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